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Abstract 

We describe a new method to formulate classical Lagrangian mechanics on a finite-dimensional 
Lie group. This new approach is much more pedagogical than many previous treatments of the 
subject, and it directly introduces students to generator matrices and their usefulness in many 
manipulations. The example of rigid body rotation, i.e. motion on the Lie group SO (3), is used 
as an example, and it is shown how to derive Euler's equations directly from the principle of least 
action. The techniques covered in this paper generalize to other Lie groups in a straightforward 
manner, which is discussed. We briefly discuss the Hamiltonian formulation of the problem as well. 
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1 Introduction 



In this paper, we derive the equations of motion for a classical system whose configuration space is a 
Lie group, using the Lagrangian and Hamiltonian formalisms. These are perhaps the most interesting 
examples of classical dynamics on non-trivial manifolds. Lie groups are manifolds which, despite 
behavior which is often topologically nontrivial at a global level, have quite simple properties locally, 
and can be completely described by a simple set of commutation relations. They appear repeatedly 
in diverse branches of physics: for example, the Lie group U(l) x SU(2) x SU(3) is the gauge group 
of the Standard Model [8]. In classical mechanics, the special orthogonal groups SO(n) relate to 
re-dimensional rigid body rotation [1]. 

While the equations of motion for systems such as rigid rotators on Lie groups are derived in 
many mechanics books as [1, 6, 5] , the derivations avoid the formalism of the principle of least action 
entirely. This misses a wonderful opportunity to introduce students to Lie groups in a context which is 
familiar: one can certainly pick up rigid bodies and throw them to observe dynamics. Furthermore, the 
derivations in traditional textbooks usually do not emphasize that concepts such as angular velocities 
and momenta take on much more subtle meanings in higher dimensions; thinking of them as vectors 
on the same footing as translational velocities and positions provides students with bad habits and 
intuition for much of physics, which increasingly takes place on manifolds beyond M 3 . For example, in 
relativity, angular momentum must be described as a tensor [5, 7]. 

The treatment of the present paper provides matrix methods for using Lagrangian and Hamiltonian 
formalisms on Lie groups, which are nice examples of manifolds where a set of coordinates is not 
inherently obvious. This not only provides pleasing theoretical closure to the theory, but it also 
emphasizes the physical nature of each of the quantities involved in the calculations. It also provides 
what is most likely the simplest example of a system where the ability to re-choose local coordinates 
at each point in time is essential. Despite these things, the set-up of the problem and the general 
procedure is similar to simpler problems in Lagrangian mechanics that are thoroughly studied. 

We first give a brief introduction to why rigid body rotation is a problem of mechanics on a Lie 
group. Then, the paper explores the Lagrangian formulation of mechanics on a Lie group, followed 
by a brief foray into the Hamiltonian formulation. Throughout the paper, the specific example of 
the group SO(3) is used, as this is the most relevant historically. Appendices are provided briefly 
discussing manifolds in mechanics, Lagrangian mechanics on a manifold, and Lie groups, and are 
necessary background for understanding this work. 

The presentation is intended to provide a sketch of how the material could be presented in a 
graduate classical mechanics course, with most derivations of interest provided directly and in a unified 
framework. The reader not familiar with manifolds or Lie groups, including the adjoint representation 
and generator matrices, is encouraged to first read the appendices, or consult a book on the subject: 
e.g., the discussion in [8]. 

2 Rigid Body Rotation 

Here we present a proof of what is frequently called Chasles' Theorem, which states that the config- 
uration space of a n-dimensional rigid body is M. n x SO(n): physically, this means translation and 
rotation, respectively. There exist many proofs of this theorem, particularly for n = 3, although most 
rely on geometric considerations [1]. Following the theme of this paper, here we use a very direct, 
computational argument. While the paper is more general than simply a discussion of rigid body 
rotation, as both the most practical and famous example of motion on a Lie group, it is appropriate 
to begin the discussion here. 
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In general, the displacement of a solid is described by u(x,t), with u 6 W 1 , x £ R ra being positions 
and t € K being time. However, for a rigid solid, the constraints on u(x, t) are strict: all lengths in the 
solid must be preserved with time, meaning that for any ti,t 2 £l and Xi,x 2 € M n , 

\u(xi,ti) - u(x 2 ,h)\ = \u(x 1 ,t 2 ) - u(x 2 ,t 2 )\. (1) 

One way to keep this constraint is to simply translate the solid: i.e., u(x,t 2 ) = u(x,ti) + a. So pick 
some point in the solid (call it x = 0) to "follow" in time, and measure all displacements relative to 
that point: i.e., u(0, t) = 0. Furthermore, because as of now we are free to label the points in the solid, 
let us define the points so that u(x, 0) = x, and thus it characterizes displacements from the initial 
position. 

Then we are left with the constraint that 

u T (x, t)u(x, t) = x T x. (2) 

We can write 

u(x,t) = A(x,t)x (3) 

where A G R nxn is a matrix. This implies that 

x T (A T (x, t)A(x, t) - l)x = 0. (4) 

A(x,t) must be x-independent: if A had an x-dependence, then we could take two ^-derivatives of (4), 
and the left term in parentheses would be x-dependent while the right term would be constant, thus 
violating the equality that the two must equal. Note that because A is the same for all points, there 
are only a finite number of degrees of freedom. 
We define O(n), the orthogonal group, as 

O(n) = {B £ E nxn | B T B = 1}. (5) 

Thinking if O(n) as a manifold, we see O(n) has two disconnected parts: one with determinant 1 
and one with determinant —1. Because the dynamics of A should be continuous, and we start, by 
construction, with A(x,0) = 1, we know that A G SO(n), the determinant 1 subgroup of O(n): 

SO(n) = {B £ R nxn | B T B = 1, det B = 1}. (6) 

Because SO(n) is also a manifold, and also has the mathematical structure of a group, because it is 
closed under matrix multiplication, it is a Lie group. 

This concludes the proof of Chasles' theorem: the configuration space consists of both the transla- 
tional part and the subsequent rotation part, and is thus M. n x SO(n). 

For the remainder of the paper, we ignore the translation on W 1 , as this is a frequently-solved 
problem. We focus on the dynamics on the Lie group, which is the topic of the present paper. 

3 Lagrangian Formalism 

We now describe how to derive equations of motion on a generic Lie group directly from the principle 
of least action. While certainly there are many different possible equations of motion on the same 
Lie group, the purpose of this discussion however is to introduce a new computational methodology 
that may be useful for many of them. We do this by presenting a derivation of the equations of 
motion of a free rigid rotator on SO(n) which naturally generalizes: we formally comment on why the 
generalization makes sense at the end. 

The general procedure is the standard one: constructing a Lagrangian and then finding the equa- 
tions of motion (Lagrange's equations); the only difference is that the procedure is kept coordinate- 
independent until the very end. 
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3.1 Constructing a Lagrangian 

In this section we formulate the Lagrangian of a system whose configuration space is SO(n). For 
simplicity, we assume that the system is not in an external potential, as this adds an unnecessary 
complication. 

Suppose we have a rigid body in n dimensions, rotating about the origin, with mass density p(x). 
If there are no external forces, the Lagrangian for the continuum motion of the body is given by the 
kinetic energy, 

L = J ^p(x)u T (x,t)u(x,t)d n x. (7) 

We have used here the notation from the previous section. Letting A € SO(n) be a matrix describing 
the orientation of the body, we can rearrange terms and, because 

ii(x,t) = Ax, (8) 

conclude that 

L = \tr[KA?A] (9) 

where we have defined the matrix 

(10) 

As of now, the Lagrangian is not expressed in terms of a coordinate chart, but this actually is helpful. 
For example, the Lagrangian certainly must have a global SO(ra) symmetry in the sense that if we 
simply rotated our frame of reference, we should not change the dynamics. That is manifest in the 
expression above, because a transformation A —■ BA for B a time independent SO(n) matrix certainly 
leaves L invariant. Furthermore, if we keep the Lagrangian in coordinate independent terms until we 
actually solve Lagrange's equations, we can pick local coordinates which both make the calculation 
more simple, elegant and easy to manipulate, and serve to emphasize the physics. 



3.2 Lagrange's Equations 

Now, we turn to the methods which can be used to solve Lagrange's equations. We repeatedly comment 
upon the connections with the mathematics to physical intuition. 

Firstly, we note that we expect that because of the global SO(n) symmetry, the equations of motion 
should only depend on the velocities in whatever generalized coordinates we choose. Secondly, we note 
that locally, we can write any matrix in SO(n) as exp[a a T a ]ylo f° r small a a . So, if we want the 
equations of motion at time t = to, we can pick coordinates so that Aq = A(to). This means that 
oc a (to) = 0. a a will serve as the set of generalized coordinates for the motion at to- One of the first 
major points of this discussion is therefore that we will exploit our freedom to re-define the a a at each 
point in time. Note that we have made no comment as to what the generators T a are, or what Aq is: 
we'll come back to this point shortly. Physically, think of this as the observer rotating his coordinate 
axes slightly at each time step, so that his coordinate axes are oriented in the same way relative to the 
body at each time. 

Given the methodology, let's express L in terms of these generalized coordinates. To do so, we use 
the fact that 



A(t ) = 1 (f^Ao) 



(a a T a A + O(a))\ t=t =a a T a A (11) 



t=to 

t=t 



where we have used the fact that since the a a vanish at t = to, only the first term contributes. (In 
general, we would have higher order corrections involving multiple generators.) This implies that, 
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taking the transpose of the above equation, using the antisymmetry of the generators and the cyclic 
property of the trace: 

L = -hr[A KA T T b T c ]a b a c . (12) 

Let's define the matrix 

M ab = -itrL4 ifA) T {T a , T b }] (13) 

where {T a ,T b } = T a T b + T b T a is the anticommutator. This matrix will show up frequently in the 
equations below. We can then write 

L = ^M bc a b a c . (14) 

Evidently M ab is like a mass/moment of inertia tensor. 

Obtaining Lagrange's equations must be done with care, but the procedure emphasizes many 
mathematical aspects which are irrelevant to the discussion on a trivial manifold where one coordinate 
chart covers the entire manifold. We begin by computing dL/da a . This derivative is in general non- 
zero, because the matrix M ab subtly carries an a a depedence in it through the Aq. Using the fact that 
(by the reasoning of (11)) 

dAo 
da a 

we obtain that using the cyclic property of the trace 



T a A , (15) 



f)T 1 

' -tr(A KA T [T a ,{T b ,T c }])a b a c . (16) 



Since 



we can simplify: 



da a 4 



\^p a rpbrpC^ rpb^rpa _J_ ^J 10, ^"^j^C Jdcdrpbrpd _|_ jabdrpdrpC (17) 



Pit i 

^ -tr(A KA T {T b ,T d })f acd a b a c = f adc M bd a b a c (18) 



da a 2 

Now, we compute the canonical momenta: 



P ^ — = M« b a b . (19) 



The time derivative here is deceptively simple: 

M ab a b . (20) 



dt 

The reason is subtle but extremely important. Because the momentum p a is only well-defined at a 
precise point on the manifold, and the coordinate chart is specifically changed at each point on the 
manifold, we are forced to use the coordinate chart at one point. At this point, M ab is a constant, 
dependent only on the position on SO(n). 
Putting together (20) and (??), we find 

M ab Co b = f adc M db uj b uj c (21) 

where oj a = a a . These are the generic equations of motion on a SO(n); indeed, their form is the same 
on other Lie groups, if the Lagrangian is quadratic in velocities. 
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3.3 Symmetries and Conservation Laws 

Before connecting the discussion to rigid body rotation on SO(3), let's investigate some consequences 
of Noether's theorem. Certainly the energy is conserved since dL/dt = 0, but here we are interested 
in consequences of group symmetries. Finding conserved quantities is always a good thing, because 
they allow us to reduce the number of differential equations that must be solved. 

As per usual, if there is any a a for which (for some choice of generators) dL/da a = 0, then p a is a 
conserved quantity. By looking at (18) this will happen for generic velocities when 

= f adc M db + f adb M dc . (22) 

for fixed a, and all b, c. We discuss the physical meaning of this in the next subsection. 

There is a more interesting class of symmetries, which come only from the group structure. If 
we define p = p a T a , ti(p k ) is an independent constant of motion for even k which are less than n, 
for SO(n) dynamics [3]. Let us prove this for k = 2. For any (semisimple) Lie group we can choose 
generators so that tr(T a T b ) ~ 5 ab , ti(p 2 ) = ti(T a T b )p a p b ~ p a p a . From the equations of motion, 

d(p P ) na-a /TCie -e fade ji ,rbd -b -c 



dt 



2p a p a = 2M ae a e f aac M bd a b a c . (23) 



since 

M ab = M ba and 

we can work in any basis, let us choose one in which M ab is diagonal (this can 
always be done, as one can show with linear algebra). This implies that 

d (P a P a ) = 2 f abc M aa M bb a a a b a c = 0, (24) 
dt J y 1 

where to obtain the last equality we used the fact that f abc = —f bac . We have used the Einstein 
summation convention whenever there are two or more indices of the same letter. In particular, we 
have shown that on a generic Lie group, the square of the "angular momentum vector" is a constant 
of motion. 

3.4 The Case of SO(3) 

Let us now turn to the case of SO (3). This corresponds to the motion of a free rigid body (no 
applied external forces) in n = 3 dimensions, and is the classic example of motion on a Lie group. We 
now exploit the fairly generic results from before, pointing out the physical meaning of the various 
quantities. 

First, note that the generators of SO(n) are antisymmetric matrices in M 3x3 . A reminder of why 
is included in the appendix. This means that the generators may be written in the form 

/ \ 

T 1 = ( 1 I , (25a) 

T 2 = | | , (25b) 

T 3 = I -1 I . (25c) 
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By computing commutators, we determine that 



[T a ,T b ] = e abc T c , (26) 

so we conclude that the structure constants f abc for this group are no different than the Levi-Civita 
tensor in 3 dimensions, e abc . 

Further manipulations with these generator matrices lead to 



{T 1 , T 1 } = ( -2 | , (27a) 
0-2 

10 

{T 1 , T 2 } = ( 1 | . (27b) 


and various cyclic permutations of the indices and respective rows/columns. This implies that, e.g., 

M 11 = j p(x) ((x 2 ) 2 + (x 3 ) 2 ) d 3 x, (28a) 

M 12 = - j p(x)x 1 x 2 d 3 x, (28b) 

and cyclic permutations of indices, in agreement with well-known results for the SO (3) rigid rotator. 
Indeed, M ab is the moment of inertia tensor. This derivation leads to important insights: quantities 
such as angular velocity and momentum, and tensors such as the moment of inertia, have indices in 
the adjoint representation. While for SO(3), the fundamental and adjoint representations have the 
same dimension (and are basically the same), this is not true in general. 

From here, we finally use the freedom to choose Aq. If we choose Aq to make M diagonal with 
entries M 1 , M 2 and M 3 , (21) becomes 

M 1 ^ 1 = uj 2 uj 3 {M 2 - M 3 ) (29) 

and cyclic permutations of indices, in agreement with the historical results. 

We also comment on the case when dL/da 1 = 0, on a chart where M ab is diagonal. (22) says that 

jabc M bb + f acb M cc = ^bc^bb _ ^-cc) = M 2 _ M 3 = ^ (3^ 

or that M 2 = M 3 . This is also an important historical result, although it is usually presented in the 
other way (e.g., in [5]): if two of the eigenvalues of M ab are equal, there is a conserved component of 
the angular momentum, which is in this case p 1 . 

It is worthwhile to explain to students that the cross product of two vectors v and w in M 3 can be 
thought of as follows: denote v = v a e a , w = w a e a where e a are the standard basis vectors for M 3 . If we 
instead think of V = v a T a and W = w a T a , where now these vectors are antisymmetric matrices, then 
we can think of v x w as equivalent to [V, W\. This follows directly from the commutation relations 
of the generator matrices. We see that the cross product is ultimately quite special and unique to 3 
dimensional spaces, in allowing us to map matrix commutation to a type of vector multiplication. 

Before moving on, it may be worthwhile to inform students that there is another Lie group with 
the same algebra as SO(3): SU(2): 

SU(2) = {B € C 2x2 | B*B = 1, det B = 1}. (31) 
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These groups, it turns out, are not isomorphic, meaning there is not a unique correspondence between 
an entry of SO(3) and an entry of SU(2): in SU(2), rotations by 2n degrees take one to —1, whereas 
in SO (3) these rotations take one back to the identity. This distinction is meaningful, as in quantum 
mechanics spinors (of spin ^) transform under SU(2) rotations, leading to the nontrivial property that 
linear combinations of spinors are not always left invariant under rotations by 2tt degrees. [9] 



3.5 Generalization to Other Lie Groups 

Let us briefly discuss why the form of the equations of motion found in this section are the equation of 
motion for "free particle" motion on any Lie group. The lowest order Lagrangian which is nontrivial 
will be quadratic in the velocities, as in (14). The coordinate dependence of L is embedded in M bc . 
Now, the canonical momenta and their derivatives are the same as before because L is the same. But 
we also expect dL/da a to be identical. This is because M bc is a tensor with indices in the adjoint 
representation, and we expect that each of them will transform in this representation under an adjoint 
representation "rotation" : 



SM bc = 5a a 



jabd J^dc _|_ jacd j^bd 



+ 0(5a 2 ). (32) 



This equation implies that we should expect for dM bc /da a , and thus dL/da a , to behave exactly as it 
did earlier. This argument provides a mathematical justification for the equations of motion, whereas 
the previous one provided a more physical interpreration. 

It is also worth mentioning that if M ab = M5 ab : i.e., the particle is a "point mass" of mass M 
which lives on a Lie group, then the equation of motion can also be thought of as the geodesic equation 
[4]. While a modification of this argument likely would also lead to a derivation of the equation of 
motion in the case of M ab not proportional to the identity, such an approach is far less intuitive to 
students who are just learning about Lie groups, as it requires associating the Christoffel coefficients, 
familiar from general relativity, with the structure constants. 



4 Hamiltonian Formalism 

We briefly comment on this formalism in the context of Hamiltonian mechanics. The Hamiltonian 
formalism of mechanics on manifolds is the study of symplectic geometry. Because there are many 
sources on this subject such as [1, 2], and the level of mathematics is significantly higher, we will 
mostly neglect a discussion of the Hamiltonian formalism. We will briefly describe it in the context of 
a Legendre transform from Lagrangian formalism. One benefit of the Hamiltonian formalism is that 
it is never unclear when to take a derivative of a quantity such as M ab , simply because the equations 
of motion are first order. 

Let (M~~ 1 ) ab be the matrix such that (M~ 1 ) ab M bc = 5 ac . Using the canonical momenta p a found 
earlier, a simple Legendre transform of 

H = p a a a - L (33) 

leads to 

H = ^(M~ 1 ) ab p a p b . (34) 

As in the Lagrangian case, the M _1 hides the dependence on the local coordinate chart a a . Nonetheless, 
using the identity {M~ l ) ab M bc = S ac and the effective calculation of dM ab / da c as done in (18), we 
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obtain Hamilton's equations to be 



oh 

dp" 
dH 



f abc (M- l ) cd p b p d = p a . 



(M- 1 ) ab p b = di 



a 



(35b) 



(35a) 



da a 



The latter equation is readily seen to be identical to the Euler equations found previously, as expected. 

5 Conclusion 

In summary, this paper presented a derivation of generalized Euler's equations for a Lie group. While 
this has been accomplished by authors in the past, the present work provides needed insight by pre- 
senting a computational trick, which is both efficient and pedagogical, to determine the equations of 
motion directly from the principle of least action. The chosen Lagrangian can thus be justified on 
either physical or mathematical grounds. The presentation is also much more appropriate for grad- 
uate students in physicis without graduate-level mathematics backgrounds in geometry; much of the 
treatment of mechanics on a Lie group is in the mathematical literature and is not presented in a form 
easily accessible by students of classical mechanics. 

The formalism provided in this paper is a bit more advanced than the usual discussions on rigid 
body motion: however, it also provides students with an introduction to Lie theory in the more 
familiar setting of finite-dimensional classical mechanics. A good understanding of Lie groups and 
their representations is crucial to understanding modern particle physics, and many physicists, such 
as the author, never see Lie theory until a study of Yang-Mills gauge theories in quantum field theory. 
Introducing the mathematics sooner would help to familiarize students more with the techniques in a 
field of physics where they are more comfortable. 

Of course, we also believe that this paper provides satisfactory methods for tackling mechanics on 
Lie groups, manifolds whose coordinate charts are challenging to obtain directly. In particular, the 
derivation of Euler's equations from first principles is far more satisfactory than typical derivations, and 
provides deep physical insight about how the nontrivial nature of the manifold leads to the nonlinear 
terms in the equations. 
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A Manifolds and Mechanics 

The configuration space of a system is the set of all possible states or orientations of our physical 
system, at any given instant in time, neglecting the rate of change of the system. We call that set Q; 
in classical mechanics, we always take Q to be a manifold; see [1] for details on the precise definition of 
a manifold. The exact details are unimportant: think of them as spaces on which one can do calculus 
(there is a notion of smoothness) and which locally look like Euclidean spaces M k for some integer k, 
but which globally may look completely different. 

Suppose that we have a free point particle of mass m moving along an infinite line in 1 dimension. 
In this case, Q = R, the set of real numbers. For our purposes, the manifold is trivial, because we 
only need to define one coordinate, x, to describe every point on Q. We say that we have defined a 
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coordinate chart, which provides an invertible map from (at least some subset of) Q to (a subset of) 
R fc (where k is the dimension of the manifold); in this case, k = 1. 

Let's consider another manifold: Q = S , or the circle (or 1-dimensional sphere, which we usually 
think of as embedded in a 2-dimensional space). If we just look at the circle in some tiny region, it 
still looks like a line; however, if we look at it "globally" it is distinct from the line: it closes back on 
itself, and is topologically distinct from R. The circle can be parameterized by an angular coordinate 
9, as is typically done in physics. However, you still only need to define 6 once to describe every point 
on the manifold uniquely. A common map sends S 1 to the subset [0, 2ir) of R. 

Now take Q = S 2 : the 2-dimensional sphere (usually thought of as embedded in a 3-dimensional 
space). Once again, this manifold locally appears flat: think of standing on the Earth and looking 
at its surface: it is hard to tell that it is curved. A typical coordinate chart here consists of angular 
coordinates (6,4>) (take eft as the azimuthal coordinate). This time, however, we run into problems. If 
9 = or 7r, then the chart to R 2 is no longer invertible: any choice of <f> works just as well. This is not 
always some abstract problem that is not important in practice. If one writes down the Lagrangian 
for a particle moving along the sphere, the equations of motion will fail at 9 = 0. In fact, it can be 
proven that there is no choice of coordinates that avoids this problem: at least one point will always 
be excluded in any chart from S 2 to R 2 . The obvious intuition is, therefore, to simply change the 
coordinates and re-orient them so that the particle avoids 9 = 0. In fact, this is a useful trick in 
Lagrangian and Hamiltonian mechanics: one can easily change the coordinate chart if the equations 
of motion will become problematic, and it is one we will exploit for good effect in deriving equations 
on a Lie group. 

It is sometimes the case in physics that Q is a manifold so abstract it is hard to visualize directly. 
We prove in this paper that Q = SO(n) for a n-dimensional rigid body with one point fixed. SO(n) is 
a subset of the set of all n x n real-valued matrices, R™*™ which is both a group and a manifold. Even 
for n = 3, the case of interest, SO(3) is a 3-dimensional space that has many non-trivial properties as 
a manifold, the most important of which, for our sake, is that it takes multiple coordinate charts to 
cover SO (3). 

While this paper is not intended to be a primer in topology, and very little topology is in fact 
needed to understand the results of this paper, it is important for students to understand how nontrivial 
manifolds arise as configuration spaces for systems that appear frequently in problems and in real life. 
The rigid body provides an excellent example of a practical problem whose configuration space is too 
abstract to even picture geometrically. 

B Lagrangian Mechanics on a Manifold 

For convenience here, as in the body of the paper, we assume that the Lagrangian is time-indepdendent. 
See [1] for more detail on many of the points below. 

Let Q be a smooth manifold, and denote by C°°(Q) the set of smooth functions / : R — > Q. We 
define the action S, a functional on Q, as a map S : C°°(Q) — > R with the property that trajectories 
which make S an extremum are the physical trajectories realized by the system. This principle of 
least action is useful because we also assume we can write the action as an integral in time over the 
Lagrangian L : TQ — > R: 



*i 

Here t\ and t<i are the starting and ending times, respectively, and TQ is the tangent bundle of the 
manifold Q, or the set of both points on Q and tangent vectors (velocities) at those points. 




(36) 
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The standard Lagrange's equations are recovered, and can be derived in the usual way, by choosing 
a coordinate chart locally on Q. If this local coordinate chart is given by {q 1 } (where the index i takes 
on a different value for each degree of freedom), we can express Lagrange's equations as 



d dL _ dL 
dt dq l dqi ' 



(37) 



This is a famous result, but we take special care to analyze it here. In particular, we remember that 
q E TgQ, the tangent fiber to the manifold at point q. This is a fancy word for the subset /submanifold 
of points in TQ that correspond to position q, and arbitrary velocity q. When we do variational 
calculus in the action to derive Lagrange's equation, we are actually evaluating dL/dqi at the point 
q + 5q, not q. Indeed, for a generic manifold, it makes no sense to talk about a tangent vector at 
point q\ as tangent at point q2, even if those points are very close to each other. On a Lie group this 
seemingly pedantic point has very real consequences. 

Having said all of this about manifolds and Lagrangian mechanics on them, it is worth emphasizing 
to students how little actually changes. The procedure of Lagrangian mechanics is identical; the only 
thing one should be careful of is that choosing a coordinate chart can be tricky, and that exploiting the 
freedom to re-choose coordinate charts (at each time) can greatly simplify the solution to the problem. 

C Lie Groups 

A Lie group G is a manifold endowed with a group multiplication operation which is smooth [2]. We 
will represent the elements of G in this paper by square real or complex matrices, so that the group 
operation is simply matrix multiplication. For a similarly brief presentation intended for physicists, 
see [8]. 

Because G is a manifold and a group, it contains the identity matrix as well as elements arbitrary 
close to the identity. We can write these elements g £ G as 



for infinitesimal e a ; in this equation, 1 is the identity matrix, {T a } is a set of matrices called generators 
for the group, and we employ the Einstein summation convention. Note that the choice of each T a is 
not unique; changing T a will just require changing e a to leave e a T a invariant. If dimG = n, the index 
a will have n possible values. 

As an example, the Lie group SO(n) can be defined as 



We can find the generators of SO(n) as follows: since we must maintain g T g = 1, we know that 




(38) 




(39) 



(1 + e a T a ) T (l + e a T a ) = 1 + e a (T a + (T a ) T ) + 0(e 2 ) 



= 1 



(40) 



which implies that 




(41) 




10 



For the purposes of this paper, it will be convenient to look at the Lie groups as slightly more 
abstract objects than as just matrix groups. In general, we can think of a Lie algebra for the generators, 
given by the commutators of the generators. If we define structure constants f abc so that 

jabcrpc (42) 

then we can define the Lie group by the f abc instead. The commutation relations above are also called 
the Lie algebra of the given Lie group. It can be shown that we can always choose f abc to be totally 
antisymmetric in its indices. A representation of the the Lie group refers to a set of matrices which can 
be written as exp[a a T a ] for some a a , where {T a } are a set of generators which obey the Lie algebra. In 
this paper, we will often find that the physically intuitive coordinates for describing (a subcomponent 
of) the system will be in the fundamental representation of the Lie group. This means that, in the case 
of SO(n), these objects are vectors in W 1 ; e.g., these are the coordinates of a point in a n-dimensional 
rigid body. We will primarily be focused on the adjoint representation. If d is the dimension of the 
group, the generators are defined as the d x d matrices with entries given by 

(n d] ) bc = r bc (43) 

in the adjoint representation. We will always use a set of coordinates a a in the adjoint representation 
to describe dynamics on a Lie group. 

References 

[1] V.I. Arnold. Mathematical Methods of Classical Mechanics (2 nd ed., 1989). 
[2] A. Cannas da Silva. Lectures on Symplectic Geometry (2001). 
[3] L.A. Dickey. Soliton Equations and Hamiltonian Systems (2003). 

[4] J.S. Dowker. "Quantum mechanics on group space and Huygens' principle", Annals of Physics 62 
(1971). 

[5] H. Goldstein, C. Poole, and J. Safko. Classical Mechanics (3 rd ed., 2002). 

[6] L.D. Landau and E.M. Lifshitz. Mechanics (3 rd ed., 1976). 

[7] L.D. Landau and E.M. Lifshitz. The Classical Theory of Fields (4 th ed., 1975). 

[8] M.E. Peskin and D.V. Schroeder. An Introduction to Quantum Field Theory (1995). 

[9] J.J. Sakurai and J. Napolitano. Modern Quantum Mechanics (2 nd ed., 2010). 



11 



